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Abstract 

^ i In order to get information about the polarized gluon distribution in a proton, we 

a I studied the electroproduction of tp' in polarized electron and polarized proton collisions 
in the framework of the NRQCD factorization approach. The value of the cross section 
dAa/dp"^ for color-octet mechanism is about 5 times larger than that for color-singlet 
one, and it might be another test of the color-octet model if the polarized gluon distri- 
bution Ag{x) is well known. Furthermore, we found that this reaction is quite effective 
for testing the model of gluon polarization by measuring the double spin asymmetry 
All for the initial electron and proton. Though the shape of p"^ distribution of All is 
quite similar for the production mechanism with color-octet and color-singlet, we can 
see a big difference in All among various models of the polarized gluon distribution 
function Ag{x). 
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The proton spin puzzle has been one of the most challenging topics in high 
energy hadron physics. As is well known, the proton spin is composed of the spin 
and angular momenta of quarks and gluons which constitute a proton. To understand 
the spin structure of a proton, it is necessary to study the behavior of the parton 
distribution functions which plays an important role in describing the spin structure of 
a proton. So far, a number of good parameterization models of the spin-independent 
parton distributions have been already obtained. In these years, we have also a large 
number of experimental data on the spin-dependent structure functions gi{x,Q'^) 
which lead to extensive study on the spin-dependent parton distributions. However, 
knowledge of the polarized gluon distribution is still poor, though many processes have 
been suggested so far for extracting information about it. 

On the other hand, since the advent of surprisingly big cross sections of prompt 
J / ip and ip' hadropro duct ion at large pt regions observed by the CDF collaboration , 
production mechanism of heavy quarkonium has been one of the most challenging topics 
in current particle physics with QCD. The CDF results were larger than the prediction 
by the conventional color-singlet model by more than one order of magnitude. J/ip 
photoproduction data at ep collider HERA experiment are also at variance with the 
prediction by the color-singlet model [Q. In order to solve these serious discrepancies 
between the experimental data and the theoretical prediction by the color-singlet model 
0, the NRQCD factorization formalism has been recently proposed as one of the 
most promising candidates. The NRQCD factorization approach separates the effects of 
short distance scales that are comparable to or smaller than the inverse of heavy quark 
mass from the ones of longer distance scales that cause the hadronization of heavy quark 
pair. A heavy quark and its antiquark pair is allowed to be produced not only in color- 
singlet but also in color-octet intermediate state at a short distance of the process, and 
subsequently the color-octet pair hadronizes into a final color-singlet quarkonium via 
emission or absorption of dynamical gluons, which is called the color-octet mechanism. 
A short distance coefficient can be computed using perturbative QCD, and a long 
distance parameter is described by nonperturbative NRQCD matrix elements whose 
values should be determined from experiments or lattice gauge theory. Unfortunately, 
the present uncertainties of the NRQCD matrix elements are not so small and thus 
the discussion on whether the NRQCD works well or not is quite controversial. There 
have been several discussions that to test the color-octet mechanism, it is important 
to study the heavy quarkonium production in various polarized reactions 0. 

In this work, to extract information about the polarized gluon density in a proton, 
we study the polarized electroproduction of ip' at large-p^ regions by taking account 
of color-octet contribution, 

e + p^ij' + X, (1) 

— * 

which could be observed in the forthcoming e-RHIC or TESLA-A^ experiments, where 
incident electron and proton are longitudinally polarized. Typical Feynman diagram 
of this process is illustrated in Fig. |T]. Since the process is dominated by photon-gluon 
fusion, the cross section must be sensitive to the gluon density in a proton and thus 
one can obtain good information about the polarized gluon distribution function. 
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Figure 1: Feynman diagram for large-pr V^' leptopro duct ion at the leading order. 
Incident electron and proton are polarized. 



An interesting observable is the double spin asymmetry All which is defined by 

where (3 is some kinematical variables. dAa represents spin-dependent cross section 
and is defined by 

dAa = - [c?(T++ — (i(7+_ + da — c/(t_+] , (3) 

where da^^, for instance, denotes the cross section with positive electron helicity and 
negative proton helicity. The asymmetry All is generally more sensitive to the mech- 
anism of hard processes than the cross section itself. This is because, since the asym- 
metry is normalized, the input parameters such as the quark mass, strong coupling 
constant, etc. which are not related to dynamics, are dropped out from the numerator 
and the denominator. 

The cross sections of large-p^ J/ ip production in ■yp collisions were calculated for 
the unpolarized case p|, [10| and the polarized case ||, |ll| . Here we use the Weizsacker- 



Williams approximation to evaluate the cross sections for the process (0). In this 
scheme we can write the electroproduction cross section in terms of the photoproduction 
cross section convoluted by the photon fiux factor in the electron, 

dAa{ep -^^' + X)= J dyAf^/,{y)dAa{yp ^ ^' + X), (4) 

where y is the energy fraction of electron carried by photon, defined hj y = Pp ■ 
Q-y/Pp ■ Pe with Pp, Pe, and being the momentum of proton, electron, and photon, 
respectively. In this calculation, we have neglected the contribution from the resolved 
photon process, which can be safely eliminated by introducing appropriate kinematical 
cuts 0. Af^/e{y) is the polarized photon distribution function in the electron, which 
is taken as |13| 



27r 



y ^iy^ \Qinax ^iy\ 



(5) 
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Figure 2: Feynman diagram for the leading order subprocess 7 + (? ^ (cc) + g a.t 
large-pr regions. 



where rUe is the electron mass. Q'^^.x maximum value of for photoproduction 

processes and we adopt Qf^ax = 1 GeV as in Ref |T^ . 



The process is dominated by photon-gluon fusion. We neglected here the con- 
tribution from the photon-quark fusion process, since its contribution to the cross 
sections is only a few percent ||^. Photoproduction cross section in the integrand of 
Eq. @ is calculated by using the subprocess cross section as 

dAai'yp + X) = J dxAg{x, Q'^)dAa{-fg ^^' + X), (6) 

where Ag{x,Q'^) is the polarized gluon distribution function in a proton. 

In the NRQCD factorization formalism, the subprocess cross section can be fac- 
torized into the short and long distance factors as 

dAa{^g ^i^' + X) = Y.^Cn{Ot,C'^'Lj)). (7) 

n 

The label n represents color and angular momentum configuration of intermediate cc 
pair. AC„ is a short distance coefficient for the polarized process and can be calcu- 
lated as perturbation series with a QCD coupling constant as- (Of g(^'^"'"^Lj)) is the 
vacuum expectation value of the NRQCD matrix element, whose relative importance is 
determined by the NRQCD velocity scaling rule. This long distance nonperturbative 
parameter represents the probability for cc pair being in the n configuration which 
evolves into physical state ip' , and can be extracted from experiments at present. 

At large-pr regions, the process is dominated by the subprocess 7* + (? — > {cc)+g 
as illustrated in Fig. ^ At the leading order, we can expect the contributions from the 
following 4 processes: 

^ + g^cc{^S,,l)+g, (8) 

l + g^cc{}SoA)+9. (9) 
l + g^cc{^S,,^)+g, (10) 
7 + (7^cc(3Pj,8)+^7, (11) 
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each of which has an intermediate cc pair with definite angular momentum and color. 
The first process is the conventional color-singlet process. The other three processes 
are color-octet processes induced by the NRQCD factorization formalism. 

For numerical values of the NRQCD matrix elements, we used the results of 



recent analysis by |14] and |15| 



(of (^^i)) = 6.70 X 10-^ 

{of C So)) = 0.75 X 10"' 

(of (^^i)) = 0.37 X 10-2 

{Of i'Po)) /ml = 0.01 X 10- 

{OfCPj)) = i2J + l){Ol 

The color-singlet matrix element given in Eq. (0) is related to the radial wave function 
at the origin, whose value is extracted from potential model calculations or directly 
from experiments. In fact, its value can be determined from the leptonic decay width 
of ip' with good accuracy. The values of the color-octet matrix elements were taken 
from the analysis on charmonium hadroproduction data. However, their values have 
not well been determined at present and various set of parameter values have been 
proposed |15]. The main source of ambiguity of those matrix elements comes from 



GeV^, 


(12) 


GeV^ 


(13) 


GeV^ 


(14) 


GeV^ 


(15) 


CPo)). 


(16) 



the higher order correction and uncertainties of unpolarized parton distribution func- 
tions. Especially, the value of P-wave color-octet matrix element (Og (^-Pq)) for the 
ip' production becomes negative in some regions, which is unphysical. For example, 
combining the data of Ref fl^ with the Type I set of Ref |]15| and the data of Ref |l^ 
with the Type II set of Ref ||15|, we obtain the value of (Og {^Pq)) as 



{OfCPo))/ml = -0.22 ± 0.13 X 10-2 GeV^ for Type I, (17) 
{Of{^Po))/ml = -0.08 ± 0.11 X 10-2 GeV^ for Type II. (18) 

As shown above, though the value of (Og {^Pq)) for Type I set is negative and unphys- 
ical, there can be a possibility of positive value of (Og {^Pq)) for Type II set by taking 
its observed error into consideration. Here we have adopted Eq. (^) as a possible 
parameter for Type II set. 

Setting Q2 = 4^2 ^j^^]^ charm quark mass = 1.5 GeV and the relevant 
center-of-mass energy y/s = 300 GeV, we calculated the spin-independent and spin- 
dependent cross sections and double spin asymmetry. The spin-dependent cross sec- 
tions are shown in Fig. |^ as a function of p"^ (left panel), and as a function of z (right 
panel), where z is defined as z = Pp -p^p'/Pp ■ q-y with p^/ being the momentum of outgo- 
ing ijj'. In this calculation, we have required the kinematical cut on pt as > 1 GeV 
for z distribution in order to suppress the contributions from the diffractive process and 



higher twist contribution. We used the GS set A [ITH] and GRSV "standard scenario" 



T7| parameterizations for the polarized gluon distribution function, and the GRV pa- 



rameterization []TB| for the unpolarized one. In these figures, the solid and dashed lines 
show the case of set A of GS and the "standard scenario" of GRSV, respectively. Bold 
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1 GS-A 

1+8 GS-A 

1 GRSV 

1+8 GRSV 




Figure 3: The spin-dependent differential cross section as a function of (left) and z 
(right). The bold lines represent the sum of the color-singlet and octet contributions, 
while the normal lines represent the color-singlet contribution only. The solid and 
dashed lines show the case of set A of GS [|T^ parameterization and the "standard 



scenario" of GRSV |17] parameterization, respectively. 



lines represent the sum of the color-singlet and octet contributions, while the normal 
lines represent the color-singlet contribution only. 

As shown in Fig. ^, the dominant contribution comes from the color-octet ^5*0 
state. We can see that the sum of the color-singlet and -octet contribution is larger 
than the color-singlet one alone by a factor of about 5 in the whole 'p^ regions. It is 
remarkable that the difference of the cross section due to these two production mech- 
anisms is larger than the one of the parameterization models for the polarized gluon 
distribution functions. Therefore, this polarized reaction could be another independent 
test of the NRQCD factorization approach. In addition, as shown in Fig. ^ (left panel), 
it might be effective even for testing the model of polarized gluon distributions, if we 
have precise values of the NRQCD matrix elements and can observe distribution of 
the cross sections precisely. 

The z distribution of the spin-dependent cross section is also shown in Fig. ^ 
(right panel). Here we do not see big difference between the color-singlet and -octet 
contribution in lower z regions. Taking account of the ambiguity of the polarized 
gluon density, we cannot distinguish between the color-singlet and -octet contribution 
in these regions. In the larger z regions the color-octet contribution rapidly increases 
and we could see a rather big difference for two production mechanisms in these regions. 
However, as is well known, the inelastic photoproduction in the NRQCD framework 
has an important issue on kinematical enhancement in the region z ^ 1 arising from 



higher order terms in the NRQCD velocity expansion |T^. It is considered that the soft 
gluon resummation becomes important in the larger z regions. The prediction for the 
cross sections by the NRQCD factorization approach in this regions should be modified 
with higher order corrections. The prediction for the leading order calculation might 
be not reliable for this region. Based on these considerations, the z distribution of 
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the spin-dependent cross section is not so effective for testing tfie color-octet model. 
As shown in Fig. ^ (right panel), it is also not good for testing the polarized gluon 
distribution in the proton. 

Next we move to the analysis on the spin correlation for inelastic ip' production 
in polarized ep collisions. The double spin asymmetries All at ^/s = 300 GeV are 
calculated and presented in Fig. § as a function of (left panel) and z (right panel). 
As seen in the p'^ distribution of All, upper two solid lines show the calculations 
for the case of GS set A parameterization. The lower two dashed lines represent the 
ones for the case of GRSV "standard scenario" parameterization. We found that the 
dependence of the production mechanisms on the p'^ distribution of All is quite small 
in the whole p"^ regions. Instead, the difference due to the models of the polarized 
gluon density is large. It becomes larger and larger with p"^. Therefore, we can rather 
clearly test the models of the polarized gluon distribution functions Ag{x) by measuring 
the double spin asymmetry for large px regions. In the z distribution, the difference 
due to the polarized gluon density is again larger than the one due to the production 
mechanisms, though the color-octet contribution dominates over the color-singlet one 
in the whole region of z. 

A similar analysis was done by Yuan et al. for the case of J/ip productions 0. 
They have discussed the color-octet contribution to J /ip production and insisted that 
the process is effective not only for the test of the color-octet mechanism but also 
for the measurement of the gluon polarization in a proton. However, only direct J/ip 
production is considered there, and they do not take account of the feedback processes 
like Ip' — *■ J /ip + X and x ~^ J/ip + X, which have been shown to play an important role 
in total J/ijj production from experiment. For ip' case, there is no such complication, 
almost all contribution come from direct production. The advantage of our analysis is 
that the analysis for ip' productions is much clearer than the one for J/ip. 

In summary, in order to obtain information about the polarized gluon distribution 
in a proton, the electropro duct ion of ip' was studied in polarized electron and polarized 
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proton collisions by taking account of color-octet contribution. As shown in the left 
panel in Fig. ^, precise measurement of dAa/dp^ for this reaction might give another 
test for the NRQCD factorization approach, if the polarized gluon distribution function 
is sufficiently established. Furthermore, if we have precise values of the NRQCD matrix 
elements, we could even test the model of the polarized gluon distribution lS.g{x) in a 
proton from it. On the other hand, d/S.a/dz is not so effective for testing the color- 
octet model and also ^g{x). As shown in Fig. ^ (left panel), the distribution of 
double spin asymmetry All{p'^) does not show a large difference of the production 
mechanism with color-octet and color-singlet. Rather, we see a big difference among 
the models of the polarized gluon distribution function Ag{x). The z distribution of 
All can be also a good test of Ag{x), if values of the NRQCD matrix elements are well 
determined. Therefore, the process is quite effective for extracting information about 
A^(x). 
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